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Abstract 

An attempt to construct the "ground state" vacuum initial data for the gravi- 
tational field surrounding two black holes is presented. The ground state is defined 
as the gravitational initial data minimizing the ADM mass within the class of data 
for which the masses of the holes and their distance are fixed. To parameterize 
different geometric arrangements of the two holes (and, therefore, their distance) 
we use an appropriately chosen scale factor. A method for analyzing the variations 
of the ADM mass and the masses (areas) of the horizons in terms of gravitational 
degrees of freedom is proposed. The Misner initial data are analyzed in this context: 
it is shown that they do not minimize the ADM mass. 
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1 Introduction 



Equations of motion of charged point particles in classical electrodynamics can be de- 
rived from field equations within the following approach: a generic state of the composed 
"particles + field" system is treated as a perturbation of the "ground state" of the field, 
uniquely determined by the positions of the particles, see pQ, [2]. In this classical approach 
the ground state is defined via a (non-local in time) decomposition of the actual field into 
the retarded (or advanced) field and the remaining "radiation field". Unfortunately, to 
decide what is the retarded or the advanced field, the entire trajectory of the particle must 
be known. Recently, it was shown that such a non-causal procedure may be avoided and 
the ground state may be defined via a conditional minimization of the (appropriately 
defined, "already renormalized" ) total energy of the "particles -|- field" system, with the 
positions and the charges of the particles being fixed (see e.g. [3], [1], [5]). Mathemati- 
cally, this leads to a simple (elliptic) variational problem for the behaviour of the field in 
a topologically non-trivial region of M.^ (exterior of the particles), where the charges of the 
particles provide the necessary boundary conditions. In this context, the ground state of 
the many particle system may be defined as the state of the field which contains ,, minimal 
amount of the radiation field" , under condition that the charges and the positions of the 
particles are fixed. 

There is a priori no obstruction against applying a similar idea to the problem of 
motion, and in particular to the two body problem, in General Relativity theory. Here, 
we have to replace point particles (elementary charges in electrodynamics) by black holes 
(elementary masses in General Relativity theory) and to consider perturbations of the 
field around a hypothetical "ground state" of the two black holes system. Such a ground 
state could be looked for as the state of the field, for which a conditional minimum of the 
total energy of the "black holes + gravitational radiation" system is achieved. Here, the 
total energy is equal to the ADM mass at infinity. This way we are led to a variational 
problem in a topologically non-trivial region of (exterior of the two holes), where the 
variations must respect: 1) the Gauss-Codazzi constraints for the field data, 2) the masses 
of the holes (assigned to each of the horizons surrounding every hole and playing role of 
the boundary conditions) and, finally: 3) the fixed positions (distance) of the holes. The 
first two conditions are technically complicated but the last one is even more difficult to 
handle because it is ambiguous: the distance between the two black holes (horizons?) 
may be measured in many nonequivalent ways. 

The idea of energy minimization is supported by the fact that, for a single black 
hole, the initial data for the Schwarzschild metric on the slice t = const, correspond, 
indeed, to the minimum of the energy (ADM mass), with the area of the minimal surface 
surrounding it being fixed (see e.g. [B], [7j, [H], [H])- Thus, the Schwarzschild initial data 
may be treated as a "ground state" of the one-body system. These data have vanishing 
extrinsic curvature and the induced three-metric is conformally flat. They lead to a static 
solution of Einstein equations. 

In electrodynamics, the ground state of the many body system is always time sym- 
metric. Indeed, any non-symmetricity means that the magnetic field does not vanish and 
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increases the energy of the state. A similar argument (non-vanishing external curvature 
could only increase the total energy of the state), together with the example of the one- 
body problem (Schwarzschild) leads us to the conclusion that the gravitational two-body 
ground state should be looked for among time symmetric field configurations only. Such 
a conjecture is supported by the fact that critical points for the ADM mass correspond 
to the initial data generating stationary vacuum space-times (see e.g. [TU], [H]). 

Similarly as in electrodynamics, where there is no static solution for the two particles 
system, in General Relativity there is no static solution for two black holes (see [12] )• The 
actual evolution of a hypothetical "two black holes -|- gravitational field" system state 
will surely display interaction between the ,, mechanical degrees of freedom" (described 
by positions and velocities of the holes) and the radiation degrees of freedom. We very 
much hope that the perturbation techniques based on the expansion with respect to the 
latter, in a phase-space neighbourhood of the ground state, could provide a good tool 
for the analysis of this evolution, at least in its initial phase. In the final stage of the 
evolution, treating the velocities on a perturbational level could become non-satisfactory. 
Then, a "dynamical ground state" would be necessary, where the time symmetry of the 
field configuration is dropped and where the velocities of the holes are encoded in the field 
boundary data on the horizons. This problem will be discussed in the future. 

In the present paper we assume, therefore, vanishing extrinsic curvature of the field 
initial data in question, which also implies vanishing of the linear and angular momentum. 
This assumption simplifies the analysis of the constraints, vector constraints becoming 
trivial. What remains is the analysis of the scalar (Hamiltonian) constraint. On the 
other hand, we do not restrict ourselves to conformally flat metrics. Such an assumption, 
used by many authors, (see [13], [H], [15]), freezes all gravitational radiation degrees of 
freedom which are contained in the conformal part of the three-metric and — as shown 
in Section [6] — is incompatible with the minimization of the energy. Our strategy is to go 
beyond conformal flatness but keeping the same topology of the exterior of the two black 
holes. As a conclusion, we assume that initial data corresponding to the ground state 
we are looking for are: 1) time symmetric (no extrinsic curvature), 2) asymptotically 
fiat, 3) the three-dimensional region S where the field data are living is topologically 
T, = M.^\Bi\B2 with two finite balls Bi and B2 being removed, and 4) the boundaries of 
these balls are minimal surfaces. 

Our goal is to construct the (momentarily static) two black holes ground state as the 
vacuum initial data for gravitational field on S (i.e. outside of the two minimal surfaces) 
for which the ADM mass (given by a boundary term at spatial infinity) is minimal within 
the set of three-metric tensors respecting: 1) the given masses of the holes and 2) their 
distance. Here, the following problems arise: 

• How do we define the mass of each black hole? This can be done in many ways, e.g. 
using quasi- local mass approach (see p^). Our choice is to assume that the mass 
of a black hole corresponds to the area of the minimal surface. 

• How do we define the distance between black holes? Instead of e.g. geodesic distance, 
we propose to use a certain global parameter defined at spatial infinity, which seems 
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to parameterize correctly the possible configurations of the two holes. 

Our analysis of the scalar constraints enables us to construct an integral identity which 
relates the variation of the ADM mass with the changes of the three-metric in the volume 
and its behaviour at the boundary. When applied to the one black hole system, our 
procedure is strongly related to the proof of the Penrose inequality ([IZ]). 

The integro-differential equations which result from our variational procedure are rela- 
tively complicated: there is little hope to be able to solve them analytically and to obtain 
an explicit formula for the gravitational ground state of the two black holes system. We 
hope, however, that an appropriate, numerical approximation will be useful for the de- 
scription of the complete two body problem. In this paper we test our methods, applying 
them to simplified situations and show that this way we easily reproduce classical results: 
the Minkowski and Schwarzschild initial data are stationary points of the ADM mass (in 
the latter case the mass of the horizon must be fixed a priori). As a by-product of our 
method we show that the two black hole Misner data do not describe the minimum of the 
ADM mass. 

The paper is organized as follows: in Section [2] we present our setup, proposed for 
description of the two-body data. We use a "2-1-1 foliation" which mimics the bispher- 
ical system of coordinates in the fiat space and we define unconstrained degrees of 
freedom which parameterize the admissible metrics. In the next Section we express the 
masses of the two interacting black holes and the total (ADM) mass of the "black holes 
-|- gravitational field" system in terms of the quantities defined in Section [21 In Section 
IHwe calculate variation of the ADM and the horizon masses in terms of the variation of 
the metric degrees of freedom. First, the calculation is performed for the conformally flat 
data and the result is apphed to the analysis of the Schwarzschild (Section [5]) and Misner 
(Section [H]) data. Section [7] is the main part of the paper, where the above analysis is 
generalized to a generic (not necessarily conformally flat) case. In Section |H] we argue, 
using symmetry requirements, that the ground state of the two black hole system has 
only one non-trivial degree of freedom. Then, we analyze in detail this simplified (yet 
physically well-motivated) case. Finally, equations for the ground state are derived in 
Section O Technical details have been shifted to the Appendices. 

2 Conformal decomposition of the metric 

Consider a three-metric g (a part of an initial data set for the gravitational field) on E, 
where S is a three-dimensional manifold with boundary, which we obtain removing two 
finite balls Ei and E2 from M^. The boundary of S consists of two disjoint connected 
(spherical) components Hi = dBi: (9S = H = HiU H2. We assume to be a Riemannian 
metric on S, asymptotically fiat in the following sense: 

g = ^^g, g = b + p, (1) 

where the conformal factor \I' is a positive function on S, b is a fiat three-metric and we 
impose appropriate decay conditions for the metric g, namely p = 0(l/r^), dp = 0(l/r^), 
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\l/ = 0(1) and (9\E' = 0(l/r^). We show in the sequel (cf. Section and Appendix [D|) that 
these decay conditions are compatible with the ground-state problem. 

We are going to use on E a "bi-spherical" system of coordinates. More precisely, we 
parameterize the one-point compactification of S by the cylinder 

S U {00} -S^xI:=W, 

where S"^ is the topological two-sphere and / is an interval [—a, b] C M. 

On W we introduce coordinates x\ i = 1,2, 3, adapted to the foliation, i.e. /i = 
G [—a, 6], and x^, A = 1,2, are spherical coordinates {ri,ip) on each of the spheres 

fi = const. We assume that Hi {H2, respectively) correspond to the value /x = —a (/i = b, 

respectively) and that infinity corresponds to the north pole [rj = 0) on the sphere /i = 0. 
Such parameterizations are subject to a three-parameter (per point) diffeomorphism 

group. To choose a specific one amongst them, three gauge conditions have to be imposed. 

To fix coordinates [r], ip) on each sphere {fi = const.}, we assume that the two-dimensional 

part of the metric g (and, therefore, also of the metric g) is proportional to the standard, 

unit, round metric on S"^: 

aAB dx^ dx^ = dr]^ + sin^ r] d(p'^ . (2) 

The above (two per point) gauge conditions give rise to the further conformal decompo- 
sition of the metric g: 

g = h'^w , 

where w is a metric on the cylinder W := S"^ x I, such that its restriction to every leave 
coincides with wab = (^ab- The entire information about such a metric is contained, 
therefore, in the following vector-density: 

■= ^detgAB g^' = aw^', (3) 

where w^^ and g^^ are the inverse (contravariant) metric tensors and a = V det aAB = sinrj. 
The remaining gauge condition, fixing the coordinate /z (and, therefore, leaves of the 
foliation {/i = const.}), may be expressed in terms of a differential equation imposed on 
the "electric field" D. In previous works [6], [7], [18], [19], conditions such as dkD'' = 
were successfully used in topologically different arrangements. Most results of this paper 
do not rely on a specific choice of the gauge. As will be seen in Section [HI some results 
simplify considerably if we assume vanishing of the mean external curvature of the leaves 
{/i = const.}, calculated with respect to the metric w (i.e. condition = 0). The 
applicability of this condition in a generic situation needs further investigation. 

Having at our disposal the product of two conformal factors: and h, we are free to 
fix arbitrarily one of them. We choose the standard value: 



which relates the flat, Euchdean three- metric b on M.^, to the cyhndrical metric w = 
d/i^ + dr/^ + sin^ r] dip"^ on W: 

b = dx^ + dy'^ + dz^ = rrrfd/i^ + dr]^ + sin^ r/dv^^l = h'^w , (4) 

(cosh /i — cos rj)'^ 

via the bispherical coordinates (/i, r], ip): 

X = cos(p — , (5) 

cosh fi — cos rj 

sinr/ 

y = sin(/? — , (6) 

cosh jj, — cos rj 

sinh/i 



cosh yU — cos 7] 

The quantity D takes in this case the following value: 

b^ = a, D^ = 0. (8) 

We shall use in parallel both of the two conformal decompositions of the metric g 
defined above. Putting 

$ := vl/v//i, 

we have: 

g = ^^w = \l/^g , 

and the following (rescaling) transformation law for the Ricci scalar holds: 

<^'n{g) = niw)-^A^<i>, 

where by we denote Laplace-Beltrami operator associated with the metric w. The 
above formula implies the following equations satisfied by the conformal factors $ and \E': 

Vdetw<i> (^A^ - $ = -^^/d^'JZ{g) , (9) 

V^vl; (^Ag - vl; = -l./d^nig) ■ (10) 

For the special case w = w and $ = ^/h = 1/ V cosh jj, — cos rj we have: g = h, TZ{w) = 2, 
TZig) = and equation Qj takes the form: 

(^A^-^Vh = 0. (11) 

However, because of singularity of the function \^ at the point /i = 0, 77 = (i.e. at 
infinity), this condition is fulfilled outside of this point, only. It is easy to check that, 
globally, the following distributional equation is satisfied, instead of fllip : 



(12) 



This means that for any (smooth) test function / G C^{W) we have: 

$ (^A^f - ^/^ = -47r/(/i = 0,r] = 0). (13) 

The Dirac delta distribution here reflects the fact that a small sphere surrounding /i = 0, 
1] = maps to a large sphere in the metric b. In the limit, point ^ = 0, t] = corresponds 
to the sphere at infinity 5*00 • 

The two ends of this cylinder, obtained for /i — > ±oo, correspond to the points 
{x,y,z) = (0,0, ±1) in the fiat, Euclidean three-space. We can interpret these points 
as two test bodies in the fiat-metric limit. The distance between these two points is stan- 
dard and equals 2. To allow arbitrary distances, there are two possible methods: 1) to 
change transformation laws ([5])-([7j), or 2) to allow an extra multiplicative factor d G M+ 
in formula (jlj). In this paper we are going to use always the last option, which implies also 
an extra multiplicative factor y/d at the right-hand side of equation ffT^ . Then g = d - h 
and the distance between endpoints /i = oo and /i = —oo equals 2d. The parameter d 
may be thought of as a "scale factor". On the other hand, the position of the interacting 
heavy bodies (black holes) in the coordinate space will always be standardized. 

In a generic (not necessarily fiat) case, the scalar constraint TZ{g) = 0, together with 
geometric identities and (fTUj) . imply the following equations: 

Ag-^^^ = 0, (14) 



- ) $ = -47rv^ (5o , (15) 



where, again, we have admitted an arbitrary scale factor \/d = \E'|^ on the right-hand 
side. (Note that w — > i& for /i, ^ 0. Hence, we may use either Vdetw or V det w in the 
integral ( JT3l) .) In the sequel, we are going to interpret 2d as a quantity parameterizing 
the distance between two black holes, even beyond the limit of point-like bodies in fiat 
space. It was already proved [20j that in the limit d —>■ oo of the Misner data, the value 
2d gives, indeed, the exact distance between the bodies. 

For our purposes, however, it is sufficient to treat d as a scale factor: its rescaling by 
a multiplicative factor implies rescaling of all the distances appearing in our problem. In 
any specific geometric situation the value of d is defined uniquely once a specific gauge 
condition for the coordinate /i is chosen. 

In a generic case of two black holes we constrain /i to a finite interval [—a, 6]. Given 
"independent degrees of freedom" of the gravitational field, described by (i.e. by the 
metric w), we want to retrieve the remaining information about the physical metric g 
from equation (fT5l) on W. The Dirac delta and the scale factor d in (fT5l) play role of 
"boundary conditions at infinity". The remaining boundary conditions on the spheres 
Hi := {/i = —a} and H2 := {/i = b} are implied by the fact that we want these two- 
surfaces to be horizons. We assume, therefore, that they are minimal in the ambient 
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three-metric g, i.e. that their mean extrinsic curvature k vanishes: = 0. We use the 
following law of transformation of k under conformal rescaling: 



= - 4Vi$ , (16) 

where k^ is the curvature calculated with respect to metric w, and is a projection 
(with respect to w) of the gradient on the unit (again in metric w), outward normal vector. 
This yields 

k\^ = ^ Vi;$L = ^ 



(17) 

In particular, if kyj = 0, we get the Neumann boundary conditions for $; for a particular, 
fiat case w = w, we have 



= 0. 



For the conformal factor \1/ these conditions may be written as follows: (compare ^Tl 
equation (36)]): 



^\oo = Vd, (^|-VI/ + ^|^(l0g V^)^ 



= 0. (19) 

H,, 



As an example of the two-body data fitting into the above frameworli: we may take the 
Misner metric [HI or its generalization see [20j for an explicit formula. The metric 
is of the form: 

and there are two minimal surfaces surrounding the singularities of the conformal factor. 
The minimal surfaces are metric spheres in the conformal metric b, hence it is natural to 
use the 2-1-1 foliation given by bispherical coordinates; the minimal surfaces correspond 
to n = ±/io (or, to fi = —a, fi = b in case of non-equal masses). The conformal metric b 
is then rewritten in terms of the cylindric metric (jlj). On the other hand, the Brill and 
Lindquist data, [15], are more difficult to handle within this framework. This is due to 
the fact that the minimal surfaces are not round spheres in the conformal metric and, 
whence, we can not use the standard bispherical coordinates to construct the foliation. 
However, it is always possible to construct a topologically bispherical foliation for these 
data, for which the metric w is no longer equal to w. 

3 ADM mass and masses of the horizons 

For time-symmetric initial data the Hamiltonian constraints reduces to TZig) = 0. Equa- 
tion ( JTOl) can, therefore, be written as: 



S^det gVg(^Vg^) = -Sv^detglVg^P - ^V^et g7^(g) , (20) 
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where |Vg\l/| is the length of the gradient of \E' with respect to g. Let us integrate formula 
fl^nj) over the whole space S. The complete divergence on the left-hand side of equation 
(!20l) yields a boundary integral. Hence we get: 



1 — >oo 



lim / ^Vg^wdetgl^^ 



8 f W^^^detgl^ 



(21) 

(8|Vg^|2 + ^27^(g)) v/d^, 

where by Vi: we denote the projection (with respect to g) of the gradient on the unit, 
outward, normal vector (analogous to V^); 5*^ are spheres of radius r (in metric b). By 
g|g we denote the pullback of the metric to S, i.e. the induced metric. We are going to 
show in the sequel that the first of the integrals yields the ADM mass: 



levrmADAf = -8 lim / W^^,/detgL. (22 



For this purpose consider a three-metric g of ADM mass m and let it satisfy the decay 
conditions we have imposed. Then, for d = 1, '^\^ = Vd = 1, the metric g has the 
following form in terms of the asymptotic spherical coordinates (r, 6, ip): 

g= (l + ^ + 0{^)y {dr^ + r^d^' + r^sin^^d^^ + Oi^)) . (23) 



Using the above form of the metric, formula fl22|) may be checked by inspection. Rescaling 
of r by an arbitrary factor d, i.e. using the Ansatz r = rd, leads to an arbitrary value 
^loo ~ following, asymptotic form of the metric tensor: 



g=(Vd+^ + 0{j,))\d7'+^d^'+^sm'^dv' + 0{j,)) . 
V 2\/dT ) 



(24) 



Simple calculations show that equation fl22l) still holds. 

Using the relation between $ and we can express the ADM mass in terms of $: 



IGvrmADAf = —8 lim 



detw|^ , (25) 



where we subtract the renormalization term ^/hV^/h (corresponding to the flat metric) 
from the term $V$. We integrate over surfaces e = const., where = /i^ + ?7^. 

We define the mass of a black hole in terms of the area of the minimal surface sur- 
rounding it0: 



^The energy niHi plays role of the lower bound in Penrose inequality (cf. [T7], [53]) and never decreases 
according to the second law of black hole physics (see e.g. [Hj or ^25j). 
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where by A we denote the two-dimensional volume element on the leaves {/i = const.}: 
A := A/det Qab- 

Given a solution $ of equation ( fTSl) with boundary conditions ( ITSl) (or, equivalently, 
a solution ^ of equation with boundary conditions ( fT9l) ) the areas of the minimal 
surfaces (i.e. the masses of black holes) are given. We have, therefore, an indirect control 
over these masses by an appropriate choice of a and h. 



4 Variations of conformally flat initial data 

We are going to prove in the sequel that the ground state of two interacting black holes 
cannot be described by conformally flat data. The search of an appropriate ground state 
must go, therefore, beyond conformal flatness. For this purpose consider a perturbation 
of the conformally fiat metric. Denote 



and decompose the degrees of freedom of the metric w as deformations of the degrees of 
freedom of w: := + 5D^, where D is given by equation ([8]). Using formula 



we have the following linearization of the above quantity: 



and, consequently, the following equation for 5$: 



A,?, - - U<l> 



1 15Z^\ ^ 

-5n^ - - — A,?,$ 



2 a 



1 



(27) 



(2J 



(29) 



We denote the right-hand side of this equation by j[w]: 



1 o 16D^ o 1 

j[w] = -57^(uO$ - A^<l> - -di 

o Z cr cr 



(30) 



Moreover, we keep the scale factor unchanged {6d = 0) and, whence, the right-hand side 
of the linearized equation ffT^ vanishes. This is why the right-hand side of fl2U|) contains 
no Dirac delta term. The boundary condition for (5$ are such that Hi remain minimal 
surfaces: 

_ d{5^) 



0. 



(31) 
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The formula fl22|) for the ADM mass can also be linearized around the metric g = h'^w: 



IQTidniADM = -8 Jii^ j (^V4;5^ + 5^V4;^) ^detgl^^. (32) 

Because the contribution of the second term under the integral vanishes, the formula may 
be rewritten as follows: 

IGnSruADM = -8}^ J (^V4<5^ - 5^V^^) ^detgl^^. (33) 

In fact, we have 5\E' — > at infinity (because the scale factor is kept unchanged: 6d = 0) 
and = O(^), which proves flHHl) . Rewriting it in terms of 5$ we obtain: 

IGTrSmADM = -8 hm j ($ V4;5$ - 6<l>Vi^^ ^det^l^^ . (34) 

We are going to express the above variation of the ADM mass in terms of variations 
SD'' of gravitational degrees of freedom. Observe that, due to elliptic equation (1^ . the 
variation (5$ of the conformal factor depends non-locally upon variations 6D^. To handle 
this non-local dependence, it is useful to rewrite the surface integral in term of a 
volume integral. Next, we shall transform the expression in such way that the dependence 
upon 6D^ becomes explicit. For this purpose we rewrite equations (J29l) and iQ in the 
following form: 

Vdet (^A^ - 1^ 5$ = Vdet w^j[w] , (35) 

Vdet liS^ (^A^ - 1^ $ = . (36) 

Subtracting these equations we get: 

Vdet w^j[w] = Vdetw (^^A^S<^ - 5$A^I>j = (^<i>V^(5$ - 5$V^<I>) . (37) 

If we multiply the above divergence by —8 and integrate it over W* := W\{{fi = 0,1] = 0)} 
then the boundary term at infinity reproduces precisely formula Finally, we obtain: 



WnSniADM = 8 f (^^V^S^ - (5$V4$) ^detw| 



H 



j $ (^A4, - ^ 5$ Vdet w . (38) 



The first integral in flHSl) vanishes because of the boundary conditions: ffT^ for $ and 
( I3T]) for 5$. Consequently, we have: 

16n6mADM = -8 i>j[w]Vdetw . (39) 
Jw* 
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In the next step we rewrite variations 5TZ and dw"^^ in j[w] (see equation flHUj) ) in terms of 
5D^. (Detailed calculations for 5TZ have been shifted to Appendix This way integral 
( !39l) may be rewritten as an expression containing variations 5D^ and their derivatives. In 
the last step we eliminate the latter using integration by parts (boundary integral vanishes 
because we assume that 5D^ vanish at the boundary). It is convenient to formulate the 
final result in terms of the following covector- valued, symmetric, bilinear form Bk{f,g): 

BA{f,9) = -\{f9)M + 2/(,A^,3) , (40) 
Bs{f,9) = I (a-i) ifg) - lcr''''lA9,B + ll39,3 , (41) 

2 

where by A we denote the Laplace-Beltrami operator on the unit sphere (in metric a). 
Then the following formula holds for an arbitrary function /: 

[ fj[w]V^= [ B,if,^)6DK (42) 
Jw Jw* 

In particular, putting / = <l>, we may rewrite formula (IH^ as follows: 



MruADM = -S 5fc($,$)(5D\ (43) 
Jw 

We want to restrict the above variation to the class of those SD^^s, for which the 
masses rriH-, i = 1,2, of both our black holes remain unchanged. For this purpose we 
are going to use the Lagrange multipliers method as soon as we are able to express the 
variations 5mj^. in terms of appropriate volume integrals. We begin with the following 
formula: 

Denote by niH^ the (unique) function satisfying equation 



- - ) m^, = (45) 



and boundary conditions 



where 5ij is the Kronecker symbol. Due to condition: V^(5$|//. = 0, integration by parts 
leads to the formula: 



IQ-KdruH^ = -8 1 {mH,Vi6<l> - ^^V^mnJ o 



(47) 
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the last equality being true because the boundary term vanishes at infinity. Hence, both 
quantities: SruADM (cf- fHHj) ) and SniHi, may be expressed by similar volume integrals of 
the form / (A^ — |) 5$, where / is a solution to, respectively: 

(^^-]]f=l' r ^"""^ (48) 

V 4y |^-47r(5o ■ Vc? for m^DM 

with appropriately chosen boundary conditions. Such / will be called respectively the 
ADM mass increase factor or the horizon mass increase factor. We already know that 
the conformal factor is the ADM mass increase factor: ffiADM = ^- Hence, we have 

IQ-ndrriADM = ~^ j ^adm - 5$Vdetw, (49) 

IQuSmu, = j {^^ - ^ S^^deiw. (50) 

A field configuration minimizing the ADM mass within the class of data with fixed 
masses of black holes must, therefore, annihilate the following form: 

lQTl{5mADM - ^^^^H) = S £ (^^(l), $) " Z/,5fc(mj,^ , $)) ^D'^ 

i 

u.mHX^] 6D\ (51) 




where Ui are Lagrange multipliers. Observe that, moreover, the above variation procedure 
respects the scale factor d which remains unchanged. We conclude that vanishing of the 
right-hand side for an arbitrary variation 6D^ is necessary if our field configuration has to 
realize the conditional minimum of the total energy of the "two black holes + gravitational 
field" system, where the masses of the holes and their distance are fixeci. 



5 Stability of the Schwarzschild initial data 

To test our method we will apply the above formulae to the fiat metric b and the 
Schwarzschild metric Qs- The fiat metric 

b = d^h^w 

describes the field configuration surrounding two "zero-mass black holes", i.e. two arbi- 
trarily chosen points of the fiat space. We have: Bk{^/h, y/h) = 0. This corresponds in 
our approach to the (weak) stability of the Minkowski initial data |26j . 

^In fact, only two variations among the three SD'^^s represent the change of the field configuration, 
whereas the third one represents variation of a gauge condition - at this point not yet fixed. 
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The Schwarzschild metric in the bispherical setting may be written as follows (cf. for- 
mula (121])): 

m y/h 



w 



Here m is the mass of the unique real horizon, whereas 2d parameterizes its "distance" 
from an arbitrarily chosen fictitious "zero-mass black hole". The radius f is a function 
of {fi,ri,(f) and parameters m and d (see Appendix [Bl for details). In this case we have: 
-Bfc($s,$s) 7^ 0, because a generic deformation 6D'' changes the mass of the hole. To 
calculate the right-hand side of (15 ip we use an explicit expression for the ttih factor: 

^ m \fh 

which may be easily verified. On the other hand, we have: 



2v^ r ' 
and, consequently: 

= v^v^+/, mu = 2f. 
The form Bk{f,g) is bilinear and symmetric, hence we get 

5fc($„ $,) - vBkifriH, $s) = 

d ■ BkiVh, Vh) + {1- 2z/)Sfc(/, /) + (2 - 2u)VdBk{Vh, /) , 

where, similarly as in the flat case, we have Bk{Vh, Vh) = 0. It is also easy to check that 
Bk{f 1 f) = 0. Hence, the above quantity vanishes for u = 1. The interpretation of this re- 
sult is following: on a class of D'^^s with appropriate asymptotics, a map D'' i— mADM[D''] 
may be defined. Constraining the map to those D'^, for which niH remains unchanged, we 
see that the linear part of rriAUMiD^] remains unchanged, i.e. the Schwarzschild metric is 
a stationary point of this map. In other words, a small deformation of the Schwarzschild 
metric which does not change the area of the horizon will not change the total ADM mass 
(cf. [6], [27!). 



6 Application to Misner data 

The two-body data, which can be easily analyzed using our method, have been proposed 
by Misner [H]. The Misner metric Qm = ^m"^ is given by a conformal factor 

^■■. = E , „, ,f , (52) 

A/cosh(/i + 2nfiQ) - COST] 



14 



defined on 5*^ x [— /io,/io]- The spheres fi = ±/io are minimal surfaces of equal mass 
m = m{fiQ,d). The above definition may be generalized to the case of non-equal masses 
[22], the explicit formulae for $ are given in pDJ. According to (E]) we have = a, 
= 0. 

Our method of variation allows us to prove that, using a small perturbation of the 
metric, one can decrease the ADM mass of the above field configurations, without changing 
the masses of both horizons and the scale factor d. (See [2H1 [22] for previous results 
obtained using different techniques.) To show this we use the fact that only two of three 
degrees of freedom are independent: an additional gauge condition fixing the foliation 
{fi = const.} may be imposed. We choose the following gauge condition: 







In this gauge we have 



15^ 



0, 



(53) 



hence, there is a function x such that we have (e^^ is the two-dimensional Levi-Civita 
symbol): 



We have also 



^ X,B 



(54) 



B35D'^ + BaSD^ = BsSD^ 



(55) 



and 



Jw* 



B.5D^ = -8 



a 



-Bas 



AB 



X 



,B 



(56) 



where Bk stands for Bk{f,^). Now 6D^ and x describe independent degrees of freedom 
and can be chosen freely. The second term in the last integral (the response to x) vanishes 
for any / of interest here because of the axial symmetry of the unperturbed metric (both 
/ and g in fj^Uj) do not depend upon the variable ip). As = 0, the response to SD^ is 
simply B^. We want to show that it is impossible to find such Lagrange multipliers ui, 1/2 
that Bsif, $m) = Bsi^m - ^ifiiH^ - U2mH2,^m) is identically zero, i.e. that Bsi^rn, ^m), 
Bsi^friHi, $m) and B3{rhH2, ^m) are linearly independent. We are unable to derive analytic 
formulae for m^i, but it is a matter of simple calculations to prove numerically this 
independence. For this purpose we approximate $m by truncating the series defining 
its value (equation (13^ in case of equal- mass data). Then we calculate the boundary 
conditions for rhn, at Hi and H2 (yU = ±/io for equal-mass data), find friH, by solving 
equation (A— l/4)m/f. = 0, and finally calculate B^^^rn, $m) and B-^ijriHi, $m)- The result 
shows that, indeed, i?3($m, $„.), B^^friHi, ^m) and B3{mH2^ ^m) are linearly independent. 
This proves that the Misner data do not minimize the ADM mass. 
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7 Variations of generic initial data 



The analysis of the Misner data presented in the previous Section shows that to construct 
the ground state we need to relax the condition w = w. This statement is also supported 
by numerical analysis of two-body data, see [30]. In this Section we describe variations 
of generic initial data, not necessarily conformally fiat. For this purpose we consider a 
generic perturbation + 5D^ of the field data and calculate the linear part of the 

corresponding perturbation of the conformal factor: $ — $ + 5$. 

To shorten the formulae we will denote V det w by ^Jw. The equations satisfied by $ 
and 5$ on W* take the form: 

A„ - ^M") » = , (57) 



where 



j[w] = l<^6n + -^{6^/w)A^<^> - -^di U6Vww'^ + Vw6w'^)dj<!>] . (59) 
8 \/w \ w 



(Observe that this is a generic form of equation 

We now calculate the coefficients -Bfe(/, $), such that 

v/^/ f A„ - 5$ = / i?,(/, $)5D^ . (60) 

w \ ° J Jw* 

We follow the procedure described in Section HI with obvious modifications. While the 
expression for in case oi w = w were quite simple, in generic case they get rather 
complicated. For example for w = w we have 5w^^ = 0, while in generic case 

the other variations get similarly complicated. Here, we give final results for Bk (the 
proofs has been shifted to Appendix [C]) . It is interesting to compare the formulae given 
below with equations fHUj) and fBTl) which hold for w = w. 
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5, 



4 (T V a 



1 (^f^ 



AB 



a- 
w 1 



cr 



D 



D 



3^2 



8 \ a 



D 



D 



GDA + (A ^ B) 



(62) 



w 



4^3 



1 

+ 4 



cr 



(^AD 



1 /v/^ 



ID, A 



2 



cr 



cr 



J- ASi) 



(63) 



In the above formulae fc^s denotes k^AB — the extrinsic curvature tensor of the {/i = 

const.} leaves with respect to w, Tabc are Christoffel symbols for the metric cr^s, 

is a normal vector, = w^'^ /yw^, and a^o* := ^(ti'*-' — M^M^)5j- In w^^. Note also the 

relation ^/w = (T\ —t. 

V 



8 Reduction to one degree of freedom 

To analyze the expressions for Bk obtained in Section [7| we need first to make some 
simplifying assumptions. In addition to the previously chosen gauge wab = o'ab we 
impose the following gauge condition: = 0, i.e. equation fl3H|) . As shown in Section [H] 
the relevant integrand is 



Bs + Ba^] 5D'-(—BAe^^ 



J 



a 



(64) 



=:B 
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The degrees of freedom are 6D^ and x- The response to 6D^, denoted by B, is a com- 
bination of i?3 and Ba, which is simpler than its components, especially if rewritten in 



terms of kab = k^AB - -jk^aAB- 



B 



1 



w 



„ABfr. r 



3 „ 



(65) 



A(/$) 



1 
4 



cr 



D 



D 



J ,B 



3\2 



C^AD 



a 



AB 



D 



^CJDA ) + (A ^ B) 



3\2 



a 



(^AD 



I) 



f^K^'^VABD 



Because of the cylindric symmetry of the problem, it is natural to assume that its 
solution (the ground state of the two black holes system) also respects this symmetry, 
i.e. that the functions do not depend upon the coordinate ip. In this case gauge 
condition 0531) together with non- singularity of the two-vector field on S"^ implies 

= 0. Moreover, it is natural to assume the axial symmetry, i.e. symmetry with respect 
to transformation cp —if. This implies that D'p = and B^ = 0. Hence, the response 
to X is zero: 



a 



-Bas 



,B 



a 



Denote 



a 



-B^e^^ 



0. 



(66) 



The function u is now the only non-trivial degree of freedom which we take into account. 
In this simplified situation we can rewrite TZ{w) as: 



n{w) = 2 + -^dAiVwa: 



2 2 
2 Am. 

u 



The response ( l65i) to 6D^ reduces to a surprisingly simple formula: 



i?3(/,<f) 



U 



U-'^,3f,3 - ^"^^^.a/.B + ^ ( A-1 ) (/$) 



and, consequently, we have: 

167r5mADM = J 53($,$)5^' 



(67) 



(6J 



(69) 



— « / a 



«-2(<|.,3)^-$^$'^ + - A-1 ($') 



Su . 
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Equation ffT^ for $ may now be written as 

d3{u-'^,3) + - (^1- Aju = -AnVddo , (70) 

where ||a is a two-dimensional covariant derivative with respect to wab = <^ab- Equation 
f lFr|) for (5$ reduces now to 

A„ - ^1 » (71) 



1 / 2 1 1 2 

= - i (m"^$^35m)^3 - (5u A$ - $'^5ti,A + -^6u - A6u 

9 Conclusions 

At this point we are able to fulfill the main goal of our paper: to formulate necessary 
conditions for the "ground state" of gravitational field around two black holes. Using 
the above (physically well-motivated) reduction to a single degree of freedom, described 
by the function u, the condition may be formulated as follows: the conformal factor $ 
satisfies equation 

ds{u~'<i>-s) + (m$II^)||a -^<I>(^1-Aju = -47:Vddo , (72) 

onW = S"^ X [—a, b] with Neumann boundary conditions at /x = —a and fj, = b. Analogi- 
cally, friHi, ^ = 1)2, satisfy equations 



(73) 



with boundary conditions (l46l) . Let us define 



2 

f = ^-^ ^i'^H, , 
i=l 

where z^i, V2 are Lagrange multipliers. To find the ground state we look for such function 
u{fi,r]) and numbers z/i, z/2 that B^i^f,^) = 0, where i?3(/, $) is given by equation (1681) . 
This condition reads: 

«-'$,3/,3 - (t^^'^aIb + ^ (a -1^ (/$) = . (74) 

We have, therefore, four equations fl7^ - fl7^ for four functions: (u, $,mj:^.) and two La- 
grange multipliers. Once these equations are solved on W, the masses of the holes can 
be read from the conformal factor $; the distance parameter is 2d. The masses can be 
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controlled indirectly because we control the parameters a and b. The asymptotic analysis 
of equation B^lf, $) = proves that the fall-off conditions for the metric g we imposed 
at the very beginning are satisfied for the ground state (see Appendix [D] for detailed 
calculations). This proves the consistency of our approach. 

Because of the high non-linearity of the problem, there is no chance to solve it an- 
alytically. We do much hope that an appropriate numerical analysis will allow to find 
solutions, which would be a good starting point for a perturbational approach. 



A Linearization of the Ricci scalar 



To derive a linear deformation of the Ricci scalar we start from the Gauss-Codazzi identity: 

2 



n{w) = R{w) + K - KabK + 



Vdet 



w 



-AWd^{KMl, + aJ)], (75) 



where k^iAB is the extrinsic curvature tensor of the leaves {/i = const.}, is its trace, 
is the normal vector, = w'^^ / \/w^, and '■= — Ml^M^dj In w^^. We have 
wab = c^AB, hence 

R{w) = 2, 

V "/ " ~ ; u u 

'"^wAB ^ 



k'ln CT k'li 



a 



w 



33 



Inserting the formulae = + 6D'', we get, in notation of Section HJ the following 
linear approximations: 




1 n'^ / 



D 



k ^ — k 



a 



^—SD^A, 
a 



AB 



wAB'^w 



Incorporating these into the formula for TZ we finally get: 



n{w)-2^--{5D^),A3 + -dA 
a a 



aa 



AB 



a 



a 



,B 



,B 



(76) 



B Schwarzschild metric — bispherical foliation 

The metric of the t = slice for the Schwarzschild solution with the mass m is (we fix 
the value of the conformal factor at infinity to be y/d): 

2y/dr 
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where r = + + (2; — -Zq)^- We introduce the bispherical coordinates 
sin ?7 cos sin ri sin sinh /i 



a; 



cosh — cos T] 



y 



cosh yU — cos r] 



cosh /i — cos 7] 



and choose such zq that the minimal surface r = 
some a > 0. This leads to the following formulae: 



m 
2d 



coincides with ji = —a sphere for 



^ m? cosh a 
Ad? sinh a 



, 2rf 

a = arsmh — 
m 



sinh 



+ 



sin^ 7] 



cosh/x — COST] ' V 4(i2 y (cosh/x — cosry)^ 



2y/d r 



f 



Tn? [ m? 

cosh /i + COST] + 2\j 1 + -^-j^ sinh + ( 1 + -^-^ ) (cosh — cosrj) 



^1/2 



C Proof of formulae (E2I) and ( 1531 ) 



In this Appendix we use notation introduced in Section [71 The first step in derivation of 
( l62l) and (l63l) is to express the linearized term ( 159|) (denoted by j[w]) in terms of It 
is a matter of straightforward calculation to see that 



w 2 D -^ 



(77) 



2.hv 



w 



6D^ 



ds Vw 9,$ 



6D^ 

5D^ 



w 



w- 



w 



6D' 
15^ 



(7J 



The remaining term, |$57^, is more complicated. As in Appendix |Al we start from the 
Gauss-Codazzi identity fl75l) . The first term in 5TZ, namely 2k^6k^, is calculated from 



1 ID^ faD'' 
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which leads to 



6k -l^fc - ' 



a 



6D^ + D 



,k 



(79) 



To hnearize the second term, k^ABk^ , we use the following formula for k^AB- 

kwAB - 



1)3^ AB, 



where F^^ are Christoffel symbols of the metric w. The expression for k^AB in terms of 
is then 



^wAB 




2 TcAB 

a a 



D 



D 



■ctda 



+ 



D 



D 



(^DB 



and the linearization can be easily calculated. The last (third) term of TZ, the diver- 
gence, involves and a^. Once these are linearized, the remaining calculations are 
straightforward. We have 



Ml 



w 



3i 



W 



33 







1 AB / ^ 



,B a 



hence 



and 



SMI 



Sat 



a 



AB 



2a 



a 6D^ 



cr 



6D" 



,B 



We have, therefore, rewritten j[w\ in terms of 5D^. Once this is done, we may read 
Bk{f, $) from the integrand ^/wfj[w]. For example, take the term 



arising from fl78p . Integration by parts yields 



w a,$ 



w* 



a 



w* \ 



a 



a 



This gives a contribution to Bs{f,^), equal ^d^^^dsf , and a contribution to BA{f,^), 
equal ^dA^d^f. The contributions to Bk{f, $) coming from all the other terms may be 
calculated in exactly the same way. Finally, after rather tedious but simple calculations, 
we obtain formulae fl62D and fl63l). 
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D Asymptotic behaviour of solutions of equation i?3(/, <l>) = 

Equation (!74l) for B3 may be analyzed asymptotically, that means near the point {fi = 
0,1] = 0) (corresponding to spatial infinity). We rewrite equation (17^ as 



$,3/,3 - «V^^<|.,^/,5 + !^ ( 1 _1 ) (/$) 



0. 



^0) 



We decompose $, / and u as 



u = 1 + u 



where $, /, u are bounded, and use Bs{\/h,\/h) = 0. Hence, vanishing of _B3(/, $) is 
equivalent asymptotically to vanishing of the following quantity: 



u 



u 



12 



(A-l)($/)-$'^/,A 



+ \/dV^,3($,3 + ^ 



- y V^V^'"($,A + + $,3/,3 + ^V^V^(A-1)($ + /) 



U 



(82) 



4 



(i/i(2M + 



d 



+ —Vd{^ + f)AVh--{2u + u^){Vh,AVh' - Vh AVh) . 
Let us expand this expression around (/i = 0,// = 0). We introduce new coordinates 



/i 

V 



—e cos 7 , 
e sin 7 . 



The asymptotic formulae 

\fh = V2e-^ + 0{e) , 7/1,3 = cos 76"^ + 0(1) , 
\fh^^ = -v/2 sin 75-2 + 0(1) , 

AVh = iV2- 3 V2 cos^ 7)e"^ + 0{£-^) = -Vh,33 

and the assumption u = 0{e) allow us to extract the lowest power of e, namely 
equation (!82l) . The only terms with are the last two terms in equation (182!) 
line). Substituting asymptotic expansions we get 



-B. 



u 



^ v^(<l> + /) (1 - 3 cos^ 7) - 1^ cos^ 7 



+ 0{e- 



(83) 

(84) 
(85) 

(86) 

from 
(the last 

(87) 
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By assumption m is a differentiable function of fi and r], hence u = e{Ci sin 7 + C2 cos 7) + 
0{e'^), Ci, C2 being constant. The functions of 7 in flHTj) are hnearly independent and the 
necessary condition for asymptotic vanishing of B3 is 

u = 0{e^) , (88) 

= 0, r/ = 0) + 7(/i = 0, r/ = 0) = . (89) 
The condition (1881) means that the metric g is of the form 

g = ^'{h + 0{r-')) , (90) 

in accordance with the asymptotic conditions imposed in Section [21 
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